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Summary. — - We present a concise review of the recent important experimental 
developments on neutrino mixing (hints for sterile neutrinos, large #13, possible non 
maximal #23, approaching sensitivity on 8c p) and their implications on models of 
neutrino mixing. The new data disfavour many models but the surviving ones still 
span a wide range going from Anarchy (no structure, no symmetry in the lepton 
sector) to a maximum of symmetry, as for the models based on discrete non-abelian 
flavour groups that can be improved following the indications from the data. 
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1. Introduction 

On the experimental side the main recent developments on neutrino mixing [1] were 
the results on 6 13 from T2K[2], MINOS [3], DOUBLE CH00Z[4], RENO [5] and es- 
pecially DAYA-BAY [6]. The different experiments are in good agreement and the 
most precise is DAYA-BAY with the result sin 2 26> i3 = 0.0890 ± 0.0112 (equivalent to 
sin 2 0i3 = 0.023 ± 0.003 or 6 13 = (8.7 ± 0.6)°). A summary of recent global fits to the 
data on oscillation parameters is presented in Table 1 [7], [8], [9]. The combined value 
of sin 0i3 is by now about 10 a away from zero and the central value is rather large, 
close to the previous upper bound. In turn a sizable 0i3 allows to extract an estimate 
of 023 from accelerator data like T2K and MINOS. There are now solid indications of a 
deviation of 023 from the maximal value, probably in the first octant [7]. In addition, 
some tenuous hints that cos Sep < are starting to appear in the data. 

A hot issue is the possible existence of sterile neutrinos: a number of hints have been 
recently reported (for a review, see [10]). They do not make yet an evidence but certainly 
pose an experimental problem that needs clarification (see, for example, Ref. [11]). 

Cosmological data allow for one single sterile neutrino but more than one are dis- 
favoured by the stringent bounds arising form nucleosynthesis (assuming that they are 
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Quantity I Ref. [7] I Ref. [8] 



AmL (10- 5 eV 2 ) 
AmL (10- 3 eV 2 ) 



sin 2 #23 
sin 2 #13 



+0.26 

-0.22 

+0.06 



7-54_ . 22 



2-43_ . 10 

0.386^;^ 
0.0241 ± 0.025 



7.50 ±0.185 

9 /I7+0.069 
^•^'-0.067 

0.30 ±0.013 

n41 +0.037 
u - 4i -0.025 

0.023 ± 0.0023 



Table I. - Fits to neutrino oscillation data. For sin 2 #23 from Ref. [8] only the absolute minimum 
in the first octant is shown 



thermalized) [12]. Actually the recently published Planck data [13] on the cosmic mi- 
crowave background (CMB) are completely consistent with no sterile neutrinos (they 
quote N e ff = 0.30 ± 0.27). The sum of all (quasi) stable (thermalized) neutrino masses 
is also bounded by cosmology. The South Pole Telescope, also using data on CMB and on 
the baryon acoustic oscillations (BAO) has claimed the result J^ m v = 0.32±0.11±? [14], 
where the question mark refers to the dependence on the assumed priors. The Planck 
experiment, also using the WMAP 9 and BAO data, finds the limit J2 m u < 0.23 at 95% 
c.L. From a completely different domain of physics, the search for neutrino- less double 
beta decay (0v/3/3), the EXO-200 experiment [15] gives an independent input on the ab- 
solute scale of neutrino masses: depending on the assumed nuclear matrix elements, they 
find m w < 0.14 - 0.38 eV. 

The MiniBooNE experiment published [16] a combined analysis of v e appearance in 
a v n beam together with v e appearance in a v^ beam. They observe an excess of events 
from neutrinos over expected background in the low energy (below 500 MeV) region of 
the event spectrum. In the most recent data the shapes of the neutrino and anti-neutrino 
spectra appear to be consistent with each other, showing excess events below 500 MeV 
and data consistent with background in the high energy region. The allowed region from 
MiniBooNE anti-neutrino data has some overlap with the parameter region preferred 
by LSND. Recently the ICARUS experiment at Gran Sasso has published the results of 
a search for electrons produced by the CERN neutrino beam [17]. No excess over the 
background was observed. As a consequence a large portion of the region allowed by 
LSND, MiniBooNE. KARMEN... is now excluded. 

Then there are v e disappearance experiments: in particular, the reactor and the 
gallium anomalies. A reevaluation of the reactor flux [18] produced an apparent gap 
between the theoretical expectations and the data taken at small distances from reactors 
(< 100 m). A different analysis confirmed the normalization shift [19]. Similarly the 
Gallium anomaly [20] depends on the assumed cross-section which could be questioned. 

Even leaving cosmology aside, if all the indications listed above were confirmed (it 
looks unlikely), then adding one or more sterile neutrinos would probably not be enough 
to satisfactorily describe all the oscillation data (see, for example Ref. [21]). There is 
in fact a strong tension between appearance (LSND, MiniBooNE, ICARUS) and disap- 
pearance (reactors, Ga anomaly) data. Thus the situation is at present confuse but the 
experimental effort should be continued because establishing the existence of sterile neu- 
trinos would be a great discovery (an experiment to clarify the issue of sterile neutrinos 
is proposed on the CERN site [22]). In fact a sterile neutrino is an exotic particle not 
predicted by the most popular models of new physics. As only a small leakage from ac- 



NEUTRINO MIXING: THEORETICAL OVERVIEW 3 

tive to sterile neutrinos is allowed by present neutrino oscillation data (see, for example, 
rcfs. [23, 24, 21, 25] and references therein), in the following we restrict our discussion 
to 3-neutrino models. 

The rather large measured value of #13, close to the old CHOOZ bound and to the 
Cabibbo angle, and the indication that #23 is not maximal both go in the direction of 
models based on Anarchy [26, 27], i.e. the ansatz that perhaps no symmetry is needed 
in the leptonic sector, only chance (this possibility has been recently reiterated, for 
example, in Ref. [28]). Anarchy can be formulated in a SU(5) ® U(1)fn context by 
taking different Froggatt-Nielsen [29] charges only for the SU(5) tcnplets (for example 
10 ~ (a, 6,0), where a > b > is the charge of the first generation, b of the second, 
zero of the third) while no charge differences appear in the 5 (e. g. 5 ~ (0,0,0)). The 
observed fact that the up-quark mass hierarchies are more pronounced than for down- 
quark and charged leptons is in agreement with this assignment. The SU(5) generators 
act vertically inside one generation, whereas the U(1)fn charges are different horizontally 
from one generation to the other. If, for a given interaction vertex, the U(1)fn charges 
do not add to zero, the vertex is forbidden in the symmetric limit. However, the U(1)fn 
symmetry (that one can assume to be a gauge symmetry) is spontaneously broken by the 
VEVs Vf of a number of flavon fields with non- vanishing charge and GUT-scale masses. 
Then a forbidden coupling is rescued but is suppressed by powers of the small parameters 
A = Vf jM , with M a large mass, with the exponents larger for larger charge mismatch. 
Thus the charges fix the powers of A, hence the degree of suppression of all elements of 
mass matrices, while arbitrary coefficients kij of order 1 in each entry of mass matrices 
are left unspecified (so that the number of parameters exceeds the number of observable 
quantities). A random selection of these kij parameters leads to distributions of resulting 
values for the measurable quantities. For anarchy (A) the mass matrices in the leptonic 
sector are totally random, while in the presence of unequal charges different entries carry 
different powers of the order parameter and thus some hierarchies arc enforced. The 
embedding of Anarchy in the SU(5) ® U(1)fn context allows to implement a parallel 
treatment of quarks and leptons. Within this framework there are many variants of these 
models: fermion charges can all be nonncgative with only negatively charged flavons, or 
there can be fermion charges of different signs with either flavons of both charges or 
only flavons of one charge. In models with no see-saw, the 5 charges completely fix 
the hierarchies (or anarchy, if the case) in the neutrino mass matrix. If Right-Handed 
(RH) neutrinos are added, they transform as SU(5) singlets and can in principle carry 
U(1)fn charges, which also must be all equal in the anarchy case. With RH neutrinos 
the see-saw mechanism can take place and the resulting phenomenology is modified. In 
Ref. [30], given the new experimental results, we have made a reappraisal of Anarchy and 
its variants within the (SUSY) SU(5) x C/(1)fn GUT framework. Based on the most 
recent data we argue that the Anarchy ansatz is probably oversimplified and, in any 
case, not compelling. In fact, suitable differences of U(1)fn charges, if also introduced 
within pcntaplets and singlets, lead to distributions that are in better agreement with 
the data with the same number of random parameters as for Anarchy. The hierarchy of 
quark masses and mixing and of charged lepton masses in all cases impose a hierarchy 
defining parameter of the order of Ac = singes with 9c being the Cabibbo angle. The 
weak points of Anarchy (A) are that with this ansatz all mixing angles should be of the 
same order, so that the relative smallness of #13 ~ o(Ac) is not automatic. Similarly 
the smallness of r = Am^ olar / Am^ tm is not easily reproduced: with no See-Saw r is of 
o(l), while in the Sec-Saw version of Anarchy the problem is only partially alleviated by 
the spreading of the neutrino mass distributions that follows from the product of three 
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matrix factors in the See-Saw formula. An advantage is already obtained if Anarchy is 
only restricted to the 23 sector of leptons as in the A^ T model (in the notation of Ref. [30] ) . 
In this case, with or without See-Saw, #i 3 is naturally suppressed and, with a single fine 
tuning one gets both #12 large and r small (this model was also recently rediscussed in 
Ref. [31]). Actually in Ref. [30] we have shown that, in the no See-Saw case, a very good 
performance is observed in a new model, the H model, where Anarchy is also relaxed 
in the 23 sector. In the H model, by taking a relatively large order parameter, one can 
reproduce the correct size for all mixing angles and mass ratios. Alternatively, in the 
See-Saw case, we have shown that the freedom of adopting RH neutrino charges of both 
signs, as in the PA^ T model, can be used to obtain a completely natural model where all 
small quantities are suppressed by the appropriate power of A. In this model a lopsided 
Dirac mass matrix is combined with a generic Majorana matrix to produce a neutrino 
mass matrix where the 23 subdetcrminant is suppressed and thus r is naturally small 
with unsuppressed #23. In addition #12 is large while #13 is suppressed. We stress again 
that the number of random parameters is the same in all these models: one coefficient of 
o(l) for every matrix element. Moreover, with an appropriate choice of charges, it is not 
only possible to reproduce the charged fermion hierarchies and the quark mixing, but 
also the order of magnitude of all small observed parameters can be naturally guaranteed. 
In conclusion, we agree that models based on chance are still perfectly viable, but we 
consider Anarchy a particularly simple choice perhaps oversimplified and certainly not 
compelling and we have argued in favour of less chaotic solutions. 

Anarchy and its variants, all sharing the dominance of randomness in the lepton 
sector, are to be confronted with models with a richer dynamical structure, some based 
on continuous groups [32] but in particular those based on discrete flavour groups (for 
reviews, see, for example, Refs. [33, 34, 35]). After the measurement of a relatively large 
value for #13 there has been an intense work to interpret these new results along different 
approaches and ideas. Examples are suitable modifications of the minimal models [36, 37] 
(we discuss the Lin model of Ref. [37] in the following), larger symmetries that already 
at LO lead to non vanishing #13 and non maximal #23 [38] , smaller symmetries that leave 
more freedom [39] , models where the flavour group and a generalised CP transformation 
are combined in a non trivial way [40] (other approaches to discrete symmetry and CP 
violation are found in Refs. [41]). 

Among the models with a non trivial dynamical structure those based on discrete 
flavour groups were motivated by the fact that the data suggest some special mixing 
patterns as good first approximations like Tri-Bimaximal (TB) or Golden Ratio (GR) 
or Bi-Maximal (BM) mixing, for example. The corresponding mixing matrices all have 
sin 2 #23 = 1/2, sin 2 6*13 = 0, values that are good approximations to the data (although 
less so since the most recent data), and differ by the value of the solar angle sin 2 #12 • 
The observed sin 2 6*12, the best measured mixing angle, is very close, from below, to the 
so called Tri-Bimaximal (TB) value [42] of sin 2 #12 = 1/3. Alternatively, it is also very 
close, from above, to the Golden Ratio (GR) value [43] sin 2 9u = —fti = Try? ~ 0.276, 

where (f> = (1 + VE)/2 is the GR (for a different connection to the GR, see Refs. [44]). On 
a different perspective, one has also considered models with Bi-Maximal (BM) mixing, 
where at leading order (LO), before diagonalization of charged leptons, sin 2 #12 = 1/2, 
i.e. it is also maximal, and the necessary, rather large, corrective terms to #12 arise from 
the diagonalization of the charged lepton mass matrices (a list of references can be found 
in Ref. [33]). Thus, if one or the other of these coincidences is taken seriously, models 
where TB or GR or BM mixing is naturally predicted provide a good first approximation 
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(but these hints cannot all be relevant and it is well possible that none is). As the 
corresponding mixing matrices have the form of rotations with fixed special angles one 
is naturally led to discrete flavour groups. 

In the following we will mainly refer to TB or BM mixing which are the most studied 
first approximations to the data. A simplest discrete symmetry for TB mixing is A4 
while BM can be obtained from S4 . Starting with the ground breaking paper in Ref. 
[45], A4 models have been widely studied (for a recent review and a list of references, 
see Ref. [50]). At LO the typical A4 model (like, for example, the one discussed in in 
Ref. [46]) leads to exact TB mixing. In these models the starting LO approximation is 
completely fixed (no chance), but the Next to LO (NLO) corrections still introduce a 
number of undetermined parameters, although in general much less numerous than for 
U(1)fn models. These models are therefore more predictive and in each model, one 
obtains relations among the departures of the three mixing angles from the LO patterns, 
restrictions on the CP violation phase 5c p, mass sum rules among the neutrino mass 
eigenvalues, definite ranges for the neutrinolcss beta decay effective Majorana mass and 
so on. Given the set of flavour symmetries and having specified the field content, the non- 
leading corrections to TB mixing, arising from higher dimensional effective operators, can 
be evaluated in a well-defined expansion. In the absence of specific dynamical tricks, in a 
generic model all three mixing angles receive corrections of the same order of magnitude. 
Since the experimentally allowed departures of #12 from the TB value, sin 812 = 1/3, 
are small, numerically not larger than 0(X C ) where Ac = sinOc, it follows that both 
#13 and the deviation of 623 from the maximal value are also expected to be typically of 
the same general size. This generic prediction of a small #13, numerically of 0{X C ), is at 
best marginal after the recent measurement of #13. 

Of course, one can introduce some additional theoretical input to improve the value 
of #13. In the case of A4, one particularly interesting example is provided by the Lin 
model [37] (see also Ref. [36]), formulated before the recent #13 results. In the Lin model 
the A4 symmetry breaking is arranged, by suitable additional Z n parities, in a way that 
the corrections to the charged lepton and the neutrino sectors are kept separated not 
only at LO but also at NLO. As a consequence, in a natural way the contribution to 
neutrino mixing from the diagonalization of the charged lcptons can be of 0(X C ), while 
those in the neutrino sector of O(Xc)- Thus, in the Lin model the NLO corrections to 
the solar angle #12 and to the reactor angle #13 are not necessarily related. In addition, in 
the Lin model the largest corrections do not affect 6*12 and satisfy the relation sin 2 #23 = 
1/2 + l/\/2 cos 6cp I sin #13 1, with Sep being the CKM-likc CP violating phase of the 
lepton sector. Note that, for #23 in the first octant, the sign of cos Sep must be negative. 
Alternatively, one can think of models where, because of a suitable symmetry, BM mixing 
holds in the neutrino sector at LO and the corrective terms for #12, which in this case 
are required to be large, arise from the diagonalization of charged lepton masses. These 
terms from the charged lepton sector, numerically of order O(Xc), would then gcnerically 
also affect 813 and the resulting angle could well be compatible with the measured value. 
An explicit model of this type based on the group S4 has been developed in Ref. [48] 
(see also Refs. [49]). An important feature of this particular model is that only #12 
and #13 are corrected by terms of O(Xc) while #23 is unchanged at this order. This 
model is compatible with present data and clearly prefers the upper range of the present 
experimental result for 813. 

In Ref. [51] we discuss three possible classes of models: 1) typical A4 models where 
813 is generically expected to be small, of the order of the observed departures of #12 from 
the TB value, and thus with preference for the lower end of the allowed experimental 
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range. 2) special A4 models, like the Lin model, where #13 is made independent of the 
deviation of sin 2 #12 from the TB value 1/3 and can be as large as the upper end of the 
allowed experimental range. In the same paper we discuss a general characterization of 
these special A4 models. 3) Models where BM mixing holds in the neutrino sector and 
large corrections to #12 and #13 arise from the diagonalization of charged leptons (in this 
case the value of #13 is naturally close to the present experimental range but #12 is at 
risk) . 

In each of the three classes of models the dominant corrections to the LO mixing 
pattern involve a number of parameters of the same order of magnitude, £. We discuss the 
success rate corresponding to the optimal value of £ for each model, obtained by scanning 
the parameter space according to a similar procedure for all three cases. We argue that, 
while the absolute values of the success rates depend on the scanning assumptions, their 
relative values in the three classes of models, provide a reliable criterium for comparison. 
We find that, for reproducing the mixing angles, the Lin type models have the best 
performance, as expected, followed by the typical A4 models while the BM mixing models 
lead to an inferior score, as they can well reproduce the size of #13 but most often fail 
to reproduce the correct value of 612- We also discuss the conditions for cos Sep peaking 
around -1. 

In Rcf. [51] we have also studied the implications for lepton flavour- violating (LFV) 
processes of the above three classes of possibilities, assuming a supersymmetric context, 
with or without See-Saw. The present bounds on LFV reactions pose severe constraints 
on the parameter space of the models (for a recent general analysis of lepton flavour 
violating effects in the context of flavour models, see Ref. [52]). In particular, we refer 
to the recent improved MEG result [53] on the /1 — > ej branching ratio, Br(fj, — » cy) < 
5.7 x 1CP 13 at 90% C.L. and to other similar processes like r — > (e or /j,)-f. In Ref. [51] 
we have studied this issue by adopting the simple CMSSM framework. While this very 
constrained version of supcrsymmetry is rather marginal after the results of the LHC 
searches, more so given that the Higgs mass is around run — 125 GeV, we still believe 
it can be used for indicative purposes as in this case. We find that the most constrained 
versions are the models with BM mixing at LO where relatively large corrections directly 
appear in the off-diagonal terms of the charged lepton mass matrix. The A4 models turn 
out to be the best suited to satisfy the experimental bounds, as the non-diagonal charged 
lepton matrix elements needed to reproduce the mixing angles are quite smaller. An 
intermediate score is achieved by the models of the Lin type, where the main corrections 
to the mixing angles arise from the neutrino sector and the non-diagonal charged lepton 
matrix elements are smaller. As for the regions of the CMSSM parameter space that 
are indicated by our analysis, the preference is for small tan/3 and large SUSY masses 
(at least one out of tuq and rri\/2 must be above 1 TeV). Given that SUSY has not 
been found at the LHC the preference for heavy SUSY masses does not pose a problem 
to these flavour models, except that, as a consequence, it appears impossible, at least 
within the CMSSM rigid framework, to satisfy the MEG bound and simultaneously to 
reproduce the muon g — 2 [54] discrepancy [55] . 
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